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A  COMPUTER  STUDY  OF  HIGH  MAGNETIC 
REYNOLDS  NUMBER  MHD  CHANNEL  FLOW 

1.0  INTRODUCTION 

This  is  a  final  report  describing  theoretical  and 
computational  studies  of  the  fully  coupled  motions  of  a  high 
speed,  high  temperature  plasma  flow  interacting  with  an  applied 
magnetid  field  and  an  external  circuit. 

The  flows  of  interest  may  be  created  by  chemical 
detonations  and  the  ensuing  process  may  be  either  inherently 
unsteady  [1]  or  quasi-steady  [2].  In  the  unsteady  process 
(Fig.  1)  the  flow  consists  of  a  hot  plasma  ("plasmoid")  formed 
between  a  driven  ionizing  shock  wave  and  its  following  contact 
surface.  The  plasmoid  is  created  by  a  sudden  release  of  energy 
into  a  driver  section  which  is  in  contact  with  a  test  gas  in 
which  the  plasmoid  propagates.  The  conducting  plasmoid  enters 
a  region  in  which  an  externally  imposed  magnetic  field  and 
electrodes  coupled  to  an  external  circuit  exist  (Fig.  1).  The 
plasma  conducts  current  to  this  external  circuit  and  is  subject 
to  Lorentz  forces  and  Joule  heating  as  it  propagates  through 
the  magnetic  field.  If  the  explosion  drive  is  a  chemical 
source,  such  a  plasmoid  will  be  of  the  order  of  5-20  cm  in 
length  in  traversing  a  magnetic  field  region  of  the  order  of 

A 

100  cm  at  velocities  of  the  order  of  10  ms  .  The  plasmoid 
may  exist  at  pressures  up  to  lk  bar  and  energies  of  up  to 
5  eV. 

In  a  steady  flow  device,  the  physical  schematic  is  the 
same  as  in  Fig.  1  except  that  a  stagnation  chamber  exists  ahead 
of  the  MHD  channel  in  which  the  explosion  driven  gases  are 
processed.  This  configuration  yields  a  hypersonic  channel  flow 
which  is  quasi-steady  over  the  flow  time  scale  in  the  channel. 

If  <Tq  ,  PqiUq'Bq  are  the  characteristic  electrical 


conductivity,  mass  density,  velocity,  and  applied  magnetic 
field,  the  flow  may  be  specified  by  an  interaction  number  i  and 
magnetic  Reynolds  cumber  rm  (in  addition  to  the  gasdynamic  Mach 
number) : 


i  = 


Rm  *  *H)r«o0 


(1) 


For  an  interaction  region  of  length  L,  the  nondimensional 
numbers  are  defined  as 


Idx 


(2) 


When  Rffl  >>  1,  the  appropriate  measure  of  the  interaction  is  the 
parameter  S  defined  as 

s  s  (Bo/^o)/<Pouo>  (3) 

where  the  spatial  average  <>  is  over  the  electrically  conduct¬ 
ing  portion  of  the  region  L.  For  a  uniform  plasmoid  of  length 

o 

a,  these  numbers  become  I  -  <TqB  a/  PqVq>  Rm  “  <J0^u0U0a* 

S  *  I/Rjjj.  In  the  case  of  a  steady  flow  situation,  the 
characteristic  length,  a,  may  be  selected  as  the  height  of  the 
duct . 


It  should  also  be  noted  that  the  square  of  the  Alfven 
2 

Mach  number,  M^,  is  equal  to  the  reciprocal  of  the  interaction 

parameter  S  (see  Appendix  A).  Thus,  electromagnetic  effects 

can  nominally  propagate  upstream  ("upstream  wakes")  for  Rm  >>  1 

only  for  S  >  1  or  i.  <  1  (see  [3]).  By  strong  interaction  we 

mean  flows  for  which  S  >  1  and  by  high  magnetic  Reynolds  number 

we  mean  R_  >  1. 
m 

The  primary  parameters  describing  the  motion  are  thus 
the  gasdynamic  Mach  number  M,  the  interaction  parameter  S  (or 
I),  and  the  magnetic  Reynolds  number  Rm.  In  addition,  a  number 
of  secondary  parameters  are  required  to  specify  the  motion. 
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These  include  the  ratio  of  electrode  length  to  duct  height,  and 
when  current  flows  to  the  external  circuit,  the  external 
current  parameter  » /ioIc/Bo  where  Ic  is  tfae  current  Per  unit 
depth  flowing  in  the  external  circuit. 


CONTACT 
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Fig.  1.  Schematic  of  explosion-driven  plasma  flow. 
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REVIEW  OF  PROGRESS 


In  these  studies  the  authors  at  STD  Research  Corporation 
have  presented  solutions  to  a  succession  of  progressively  more 
complex  initial  value  and  boundary  value  problems  building  up 
to  the  strong-interaction,  high  magnetic  Reynolds  number 
motion.  Thus  in  [4],  the  strong  interaction,  high  magnetic 
Reynolds  number  motions  were  studied  in  the  quasi-one- 
dimensional  approximation.  In  that  work  several  important 
phenomena  were  revealed  for  the  first  time.  It  was  shown  that 
under  conditions  of  strong  interaction  and  moderate- to- high 
magnetic  Reynolds  numbers  the  interaction  zone  progressively 
decouples  from  the  shock  front  and  evolves  into  a  magneto¬ 
hydrodynamic  discontinuity  (suitably  dispersed  depending  upon 
the  magnetic  Reynolds  number)  which  moves  at  the  particle  speed 
of  the  fluid  rather  than  at  the  shock  front  speed.  It  was 
further  shown  in  (4J  that  upstream  running  shock  waves  will  in 
general  be  created  in  the  interaction  zone  under  conditions  of 
large  interaction  parameters. 

In  [5]  the  authors  considered  the  two-dimensional  motion 
under  conditions  of  low  to  high  magnetic  Reynolds  numbers  but 
under  conditions  of  weak  Interaction.  In  this  case,  the 
electrical  motion  is  that  which  results  from  a  two-dimensional, 
high-speed,  turbulent  flow  proceeding  in  a  constant  area  duct; 
since  the  interaction  parameter  is  vanishingly  small  the  fluid 
motion  is  not  Influenced  by  the  electrical  motion.  Since  there 
are  no  geometrical  nonuniformities  or  magnetohydrodynamic 
forces  operating  on  the  fluid,  there  are  no  shock  wave 
interactions  set  up.  The  principal  fluid  nonuniformities  are 
the  turbulent  boundary  layers  generated  on  the  duct  walls.  In 
these  studies  the  detailed  eddy-current  motions  in  the  non- 
uniform  regions  of  the  applied  magnetic  field  were  revealed; 
but  more  Importantly,  the  current  distribution  within  the 
conducting  fluid  resulting  from  the  electrode  system  and  the 
external  circuit  current  were  also  revealed.  In  particular  it 
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was  shown  that  for  magnetic  Reynolds  numbers  in  the  range  10  < 
R_  <  20  the  current  flowing  between  the  electrodes  is  convected 
fully  downstream  of  the  electrode  gap.  As  a  result,  the 
Lorentz  forces  which  would  act  on  the  fluid  exist  downstream  of 
(and  not  within)  the  interelectrode  gap.  It  was  further  shown 
that  even  for  electrode  systems  which  are  maintained  such  that 
the  electrode  edges  do  not  extend  into  the  fringe  field  of  the 
applied  magnetic  field,  there  is  still  considerable  coupling  of 
the  fringe  field  induced  eddy-cells  and  the  electrode  system, 
particularly  when  the  external  circuit  is  open. 

In  [5]  the  authors  studied  a  quasi- three-dimensional 
description  of  the  strong  interaction-high  magnetic  Reynolds 
number  regime.  There  it  was  shown  that  the  Lorentz  forces 
concentrated  in  the  magnetic  boundary  layers  under  strong 
interaction  conditions  would  stall  these  boundary  layers  and 
establish  a  recirculation  cell  system  in  the  duct  with  a 
central  jet  of  plasma  which  would  be  relatively  free  of  Lorentz 
forces.  Because  of  the  restrictions  of  the  quasi- three- 
dimensional  model,  the  ref lected-upstream  moving  shock  waves  as 
well  as  the  details  of  the  recirculation  cells  cannot  be 
described.  Only  the  onset  of  stall  and  recirculation  can  be 
predicted. 
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APPENDIX  A 


THE  MATHEMATICAL  DESCRIPTION  OF  MAGNETOGASDYNAMIC  FLOWS 
AT  STRONG  INTERACTION  AND  HIGH  MAGNETIC  REYNOLDS  NUMBER 


2.  THE  MATHEMATICAL  DESCRIPTION  OF  MAGNET OGASDYNAMIC 


FLOWS  AT  STRONG  INTERACTION  AND  HIGH  MAGNETIC  REYNOLDS  NUMBER 
2.  i  Fluid  Conservation  Laws 

We  may  describe  the  fluid  in  terms  of  its  mass  density,  p, 
velocity  U,  and  internal  energy  e  .  We  describe  the  electromagnetic 
effects  in  terms  of  the  electric  field  E  and  magnetic  field  S’.  These 
variables  are  considered  to  be  general  functions  of  space  He  and  time  t. 

The  conservation  laws  for  mass,  momentum,  and  energy  are 


+7  •  (pU)  =  0 


^  (PU  )  +  7  •  (pUU  )  =  7  •  n  +JXB 


^  £P(e  -  U2/2)J  +  7  •  +U2/2)  puj=V-(ir-u)-V-q+J-E 


In  the  conservation  laws,  «  is  the  total  pressure  tensor  and 
q  is  the  heat  flux  vector.  This  system  of  conservation  laws  is  completed 
in  the  limit  of  infinitely  fast  kinetics  by  the  kinetic  and  caloric  equations 
of  state 

p  =  P<P  ,  *  )  (4) 

e  =  c  (p,T)  (5) 


where  p  is  the  isotropic  part  of  the  stress  tensor  ff  and  T  is  the 
temperature.  For  a  general  fluid  the  state  equations,  Eqs.  (4),  (5)  cannot  be 
explicitly  given  but  are  embedded  in  the  general  statistical  mechanical 
description  of  the  equilibrium  thermochemistry  of  the  system. 
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In  the  case  of  a  perfect  gas  with  particular  gas  constant  R 
and  specific  heat  ratio  y  explicit  formulae  may  be  given: 


p  =  pRT 

(6) 

e  =  (y-D^RT 

(7) 

p  =  p<Y-Oe 

(8) 

2.2 

The  Electromagnetic  Contributions 

The  electrical  equations  (consisting  of  the  Maxwell  equations 
anA  the  generalized  Ohm's  law)  govern  the  electric  and  magnetic  fields 

E  ,B  and  the  conduction  current  density  J  .  In  the  hydromagnetic  limit 

these  are 

VxE  :  • 

ot 

(9) 

V  x  B  =  p-qJ 

(10) 

7  .  B  =0 

(ID 

7  =  <r  (e  +  Uxb)+7k 

(12) 

where  is  the  thermal  diffusion  flux  vector, 

“K  =  cK 


and  K  is  given  by  l 1 ] 


K 


9(1,7Te  +  0 


(2) 


7T  XB 
e 


+  9^  (VTeXB)XB 


Vpa+(3 


(2) 

a 


V  pa  X  B  +  0 


(3) 


( *7  p  X  B 
ra 


)  X  B 


(13) 
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The  9^1  •••  and  p^  •  •  •  are  transport  coefficients  defined  in 
[ljf  is  the  electron  temperature  and  the  subscript  a  denotes  a 
plasma  component. 

The  contributions  to  the  momentum  and  energy  of  the  system 
by  the  electromagnetic  field  are  contained  within  the  Lorentz  force, 

J  XB  and  the  Lorentz  power  J  •  E.  The  Lorentz  force  may  be 
represented  in  terms  of  the  Maxwell  stress  tensor  T,  as 

J  XB  =  V  •  T  (14) 

where  the  Maxwell  stress  tensor  is  defined  as 

T  =  (b  B  -  B2/2  T)  (15) 

Correspondingly  the  Lorentz  power  may  be  represented  in  terms  of  the 
Poynting  flux  S  and  the  electromagnetic  energy  density  e^ 

— 

J  •  E  =  -V  •  S  -  g-p2-  (16) 

The  Poynting  flux  S  is  defined  as 

T  =  u Q" 1  E  X  B  (17) 

while  the  electromagnetic  energy  density  is 

em  2  »o‘  b2/2  <l8> 

Let  us  expand  the  Poynting  flux  in  terms  of  J  ,  B  through  the  use  of 
the  Ohm's  law 


T  =  <t(e+uxb+k) 


We  have 


1  =  Tq  (sxb)  =  (^0<r)  J  XB  -Hq1  (u  XB)  XB  -^‘kxB 

The  term  7  X  B  is  simply  7  •  T  .  The  term  (U  XB)XB  is  readily 
shown  to  be 

(UXB)XB  =  U  •  (BB)  -U!(B  -B) 

which  can  be  rearranged  to 

(UXB)XB  =  U  •(BB  >  BZ/lT )  -  U  (B2/2) 

The  Poynting  flux  is  therefore  represented  as 

r  =  e  U  -  U  •  F  +  rj7- T-nTKXB  (19) 

/\  - 1 

where  n  s  (nQ<r)  is  the  magnetic  diffusivity.  We  note  that  the  Poynting 
flux  can  be  decomposed  into  four  constituent  parts:  (a)  a  purely  convected 
flux  of  electromagnetic  energy  carried  by  the  motion  of  the  medium 
(emU  );  (b)  a  power  flow  represented  by  work  done  per  unit  time  by  the 
Maxwell  stresses  acting  on  the  moving  medium  -ff"  -T  ;  (c)  a  diffusive 

flux  of  electromagnetic  energy  driven  by  gradients  of  the  Maxwell  stress 
tensor  (Hi7*r);  and  (d)  a  power  flow  represented  by  work  done  per  unit 
time  by  the  thermal  diffusion  flux  (-r|  x  B). 

The  electrio  and  magnetic  fields  and  currents  may  be  expressed 
in  terms  of  vector  potentials  A  .and  scalar  potential  9  a  a 


b"  =  7  X  A 

(20a) 

?  ■  -**  -  w 

(20b) 
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2.3 

Fluid- Electrical  Svstem 

viscous. 

The  mass,  momentum,  and  energy  equations  for  the 
hydromagnetic  system  may  now  be  expressed  as 

general, 

It  ‘  *  5  0 

(21) 

+  v  •  r  =  7  •  t* 

at 

(22) 

•  ~H  =  V-  (IJ  *"t)-  V-  (rjV-7)  -  7-‘q  +  V'(n\xB 

(23) 

In  the  above  M  -  pU  is  the  momentum  density  and  T  is  the  total  fluid 
and  electromagnetic  momentum  flux 

V  =  pUlT  +  pi  -  F  (24) 


The  total  energy  density  is 


e  =  p(€  +  UV2)  +  e 


m 


(25) 


and  g  is  the  total  enthalpy  flux  vector 

7  =  [<e  +  P)T  -  r]  -  u 

The  fluid  stress  tensor  V  has  been  decomposed  into  a  pressure  p  and 
viscous  stress  tensor  r  where 


(26) 


p  s  ■  j  Trace  [it  ) 


(27) 


and 


It  =  -pi  +  T 


(28) 


We  may  define  the  electromagnetic  pressure  pm  as  the  mean 
normal  compressive  Maxwell  stress: 
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(29) 


Pm  =  "  3  Trace(7  > 


Expressed  in  terms  of  the  magnetic  field  intensity  IT,  the  pressure  is 

Pms3(B2/2^o)  (30) 


or  expressed  in  terms  of  the  electromagnetic  energy  *  density  e 


Pm  ~  3  *m 


m 


(31) 


The  Maxwell  stress  tensor  may  be  decomposed  into  a  mag¬ 
netic  pressure  pm  and  a  Maxwell  stress  deviator  from  isotropy  18 


*  s  **  "  Pm1 


(32) 


Equation  (32)  may  be  thought  of  as  the  defining  equation  for  the  Maxwell 
stress  deviators  T  The  mass,  momentum,  and  energy  equations  may 
now  be  rewritten  in  alternative  form  as 

+  V  •  M  =0  (33) 

+  V  •  G  =  7-  F*  ♦  V  •  7  (34) 

||  +  7  •  H  =  V  •  (U’V  +  7  -(U-7)  -  7  -  cf  -  7  -(nv-  n  +  7  *(nTKxB)  (35) 

In  the  above,  G  is  the  total  momentum  flux  with  only  the  magnetic  pressure 
included  as  the  magnetic  contribution 

G  =  pUU  +  (p  +  pm)  T  (36) 

while  H  is  the  total  enthalpy  flux  consisting  of  both  fluid  and  electro¬ 
magnetic  pressure  and  energy  contributions. 

H  =  [e  +  (p  ♦  pm)]  U  (37) 
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In  terms  of  the  total  energy  density  e,  the  momentum  density 
14,  and  the  electromagnetic  energy  density  em,  the  state  equations  (6)  - 
(8)  become 


p  *  "  *m‘#2/2pJ  (38) 

T  =  p/pR  =  ^[e  -  «m  -  «2/2p]  (39) 


Let  us  now  consider  the  transformation  of  the  electrical 
equations  (9)-(12)  into  more  useful  forms.  Combining  Cqs.  (9)“(12)  we 
obtain  the  governing  equation  for  the  magnetic  induction  B: 

-  7  x  (U  xB)  »  -7  x  (n7  x  S')  +  V  x  K  (40) 

o  t 

We  note  that  given  the  magnetic  induction  B  ( x,  t)  governed  by  Eq.  (40) 
one  immediately  has  specified  the  Maxwell  stress  tensor  T  and  the 
electromagnetic  energy  density  «m-  Further,  the  current  density  J  is 
determined  from  B  as 


T  *  Kg1 7  X  B 


and  the  electric  field  E  as 

E  s  -if  x  B  +  1)7  X  B  •  K 


2. 4  Viscous  and  Heat  Conduction  Effects 

Let  us  now  make  some  observations  about  the  viscous  stress 
tensor  T  and  the  heat  flux  vector  "q  .  The  Navier-Stokes  moments  of 
the  Boltzmann  equation  yield  kinetic  theory  forms  of  these  quantities: 


rL  *  2p(7U)0 


(41) 


<*L 


\7T 


(42) 


V 


In  the  above,  k  and^X  are  the  coefficients  of  viscosity  and  thermal 

conductivity  and  (VU)Q  is  the  symmetrized,  traceless  velocity  gradient 
tensor: 


•—  ,  /9U.  91L,  .  3U. 

,VU)Oij  *  i  ’  3^)  -  I  6ij 


We  denote  the  stress  tensor  and  heat  flux  vector  with  a  subscript  L  to 
denote  that  these  are  laminar  quantities.  If,  on  the  other  hand,  we 
interpret  the  fluid  variables  p,U,e,T,  •••  as  turbulent  mean  quantities, 
then  T  and  'q  contain  turbulent  contributions  due  to  the  turbulent  velocity 
and  enthalpy  correlations.  Hence,  the  complete  stress  and  heat  flux 
fields  for  a  turbulent  hydromagnetic  medium  are 

7  s<po'  u'>+ ~L 

”q  =  <  pU  h  )  +  +  Tl 

where  U  ,h  are  the  turbulently  fluctuating  velocity  and  enthalpy  and  (  ) 
denotes  an  ensemble  average.  A  detailed  higher  order  closure  theory  for 
the  turbulent  contributions  (pU  U  ),  (pU  h)  is  given  by  Demetriades, 
Argyropoulos,  and  Lackner  [2],  The  radiative  heat  flux  is  q^.  Since 
the  optical  depths  in  dense,  explosion  generated  plasma  are  so  small, 
the  radiative  heat  flux  is  only  important  in  layers  near  the  plasma  surface 
of  the  order  of  the  radiation  free  path. 


Let  us  consider  the  nondimensionalization  of  the  fluid  equations 
(33)  -  (35)  and  the  electrical  equation  (4  ).  For  this  purpose  let  us 
specify  characteristic  values  of  the  variables  as  p0,Ug,«Q,pg,  TQ,  •  •  •  as 
well  as  magnetic  field  B q.  We  define  a  characteristic  length  L  and 
characteristic  time  Iq  =  L/Uq.  We  indicate  nondimensional  variables 
with  (  ). 
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The  fluid  and  electrical  conservation  laws  then  become 


^  +  V  -  M  *  0 
3t 


•&JL  +  V  .  f  = 
at 


R^7-  2MVU)0 


%+7-fi  *  R^V-^JIU  •  (VU)o]+ [(v-DM^pJ  V-(WT)  - 


-  SR^V-(nV-r)  -  2S7.(nJKxB) 


(43) 

(44) 


(45) 


^•Vx(uxB)=  "Rm{^  X  (nVxB)  -  7x  (n  J^)j 


(46) 


B  =  V  X  A 


E  *  -7*- 

at" 


In  the  form  in  which  the  Maxwell  stress  tensor  is  decomposed  the  left  hand 
sides  of  the  momentum  and  energy  equations  take  the  forms 


87.?.  .2# 

+  7  .  f 

at 

at 

^  +  7-  H-S7  *(U  •  TJ  =  2* 

+  7  •  i? 

at 

at 

The  nondimensional  variables  are  defined  as 


t  •  *  UQt/L  7  =  L7 

U  *  uT/u0 

P  *  P/P0 

"  *  */p0u0 


to 


T  =  T/Tq 

e  =  P  [yP^M2]  +  U2  j-  +  S  em 

«  =  ^/[(y-D^RTJ 

®m  =  <b2/2P0)/(B2/2?0) 

E  =  E/U0B0 

B  *  B/Bq 

A  »  A/(Bq/L) 

*  *  */(U0B0L) 

J  »T/(P0mb0/l) 

!V¥0B0> 

G  *  pU  U  +  [ftM2)1  5  +  SpJ  T 
f  a  P  U  u  +  (vM2)' 1  pT  -  sr 

Pm  =  Pm/(Bo/2Po> 
h  =  [;♦  (ym2)  ;  +  sjm]  u 

h  •  [(«  +  (ym2)1J)T  -  sr]-u 
?  «  r/(B2/2p0) 

*  *  ^Aq  h  *  ^/n0  n  *  n/n0 

where 

xo  *  ^Po'V'  »*o  s  ^(Po*To)'  '’o  ■|*<Po»To)- 

The  state  equations  (38 )f  (30)  in  nondimensional  form  are 


tl 


p 


Y(Y-l)M2je 


P 


It  can  be  seen  that  the  general  viscous  hydromagnetic  equations 
contain  six  fundamental  nondimensional  parameters.  These  are 


Mach  number 
Interaction  parameter 
Magnetic  Reynolds  number 
Viscous  Reynolds  number 
Viscous  Prandtl  number 


V<W(Ao,1/2 

(Bj/2S0)/(P0U|) 

<uox./n0) 

CoU0L/'“0 

«V<>A0> 


We  note  that  the  Alfv£n  speed  is  defined  as 

1/2 

CA  =  (b2/^op) 

and  the  Alfv£n  Mach  number  is  M^  =  U/C^.  Hence  the  interaction  param¬ 
eter  S  is  also  twice  the  reciprocal  of  the  square  of  the  Alfven  Mach  number. 

For  flows  in  the  absence  of  viscous  and  diffusion  effects 
Rm  -*  «,  R#  —  oo  we  obtain  the  inviscid  hydromagnetic  equations: 

11  +  V  •  1  =  0  *  (47) 

at 

iL£-  +  7.j*G  ■  S  f,  1  8  0  (48) 

at  J 
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^  +  V-  J  H  -  U  •  r  *  1=0  (49) 

at  l  "  J 


^--VX(UXB>  =  0  (50) 

at 


The  jump  equations  across  hydromagnetic  shocks  immediately  follow 
from  Eqs.  (47) -(50). 


(51) 


(52) 


(53) 


(54) 


where  {  }  denotes  the  difference  in  the  quantity  across  the  shock  surface 
and  IT  is  the  normal  to  the  shock  surface . 

Since  the  electrical  conductivity  achievable  in  nonideal  plasma 
is  large  but  finite,  the  magnetic  Reynolds  numbers  are  not  infinite  but 
perhaps  vary  in  the  range  1  —  Rm  —  20. 

In  this  range,  the  appropriate  system  is  the  inviscid,  finite 
conductivity  hydromagnetic  systems 

^  +  7  •  M  =0  (55) 

at 

-24-  +  7  •  f  =0  (56) 

at 
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(5?) 


^  +  V  •  5  =  -SR'1  v  •  in  7  •  ?) 
at  m 


3  R  —  —  —  ^  ~~ 

“  •  v  X  (U  X  B)  =  -R  7x(nVxB) 

at  m 


(58) 


or  in  the  alternative  form,  the  momentum  and  energy  equations  are 


+  V  •  G  =  s  7  •  T 
at  * 


(56a) 


2.6 


%  +  v- H  =  S  V-  (U  •  X  )  -  SR^V  •(  rj  V*  T) 

at  *  m 


Applied  and  Induced  Fields 


(56b) 


n(0) 


Let  us  separate  the  magnetic  field  B  into  an  applied  portion 
B  y'J’  sustained  by  currents  external  to  the  plasma  and  plasma  induced 
portion  B  ^  which  results  from  currents  flowing  within  the  plasma: 

B  =  B(0)  +  B(i) 

The  induction  equation,  Eq.  (58),  then  becomes 

T(i) 


(57) 


»  r^(i)  «.  —  -r(i) 

— - Vx(U  x 

at 


r(0) 

3  B 

where  — 3—  is  denoted 

at 


“w  -if-  —  ~  ~ ~  'i 

B  = -Rm-jVx  (rjVxB  )-Vx(r)JK)j- 

-  r(0)  t(0) 

•R ‘7x(hVxb  )  -  B 
m 


(58) 


B  (0) 
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Abstract 


The  magnetohydrodynamlc  behavior  of  shock¬ 
generated  pulses  of  plasma  ("plasmolds*)  In 
encounter  with  an  applied  magnetic  field  is 
studied.  Reflection  and  transmission  of  the 
plasmold  shock  front  is  revealed  as  well  as  the 
existence  of  Joule  heated  zones  which  are 
convectlvely  unstable.  It  is  shown  that  strong 
interaction  plasmolds  are  not  delimited  by  shock 
front  and  following  contact  surface;  rather,  they 
develop  their  own  internal,  evolving  structure 
under  the  mutual  influence  of  self  Joule  heating, 
and  Induced  magnetic  field.  As  a  result,  the 
interaction  zone  behind  an  ionizing  shock  wave  will 
progressively  decouple  from  the  shock  front  and 
follow  a  trajectory  given  by  the  particle  motion  of 
the  gas  rather  than  the  wave  motion  of  the  shock 
front . 


I.  Introduction 


In  the  present  work  we  study  the  behavior  of 
shock-generated  magnetohydrodynamic  flow  over  the 
range  of  weak  and  strong  magnetohydrodynamlc' 
interaction  and  low  to  high  magnetic  Reynolds 
numbers.  Such  a  flow  consists  of  a  hot  plasma 
■plasmold"  formed  between  a  driven  ionizing  shock 
wave  and  its  following  contact  surface.  The 
plasmold  is  created  by  a  sudden  release  of  energy 
in  a  driver  section  which  is  in  contact  with  a  test 
gas  in  which  the  plasmold  propagates.  Such  a  flow 
may  be  driven,  for  example,  by  the  use  of  focused 
chemical  explosives1  *J. 

The  conducting  plasmold  enters  a  region  in 
which  an  externally  Imposed  magnetic  field  and 
electrodes  coupled  to  an  external  circuit  exist 
(Fig.  1).  The  plasma  conducts  current  to  this 
external  circuit  and  is  subject  to  Lorentz  forces 
and  Joule  heating  as  it  propagates  through  the 
magnetic  field.  If  the  explosion  drive  is  a 
chemical  source,  such  a  plasmold  will  be  of  the 
order  of  9-20  cm  in  length  in  traversing  a  magnetic 
field  region  of  thjj  order  of  100  cm  at  velocities 
of  the  order  of  10  ns”  .  The  plasmold  may  exist 
at  preaaures  up  to  1  k  bar  and  energies  of  5  eV. 

If  ej,.  P-.U  are  the  characteristic  electrical 
conductliAty,  mass  density  and  velocity  within  the 
plasmold,  the  flow  may  be  specified  by  an  inter¬ 
action  number  per  unit  length  1  and  magnetic 
Reynolds  number  per  unit  length  r  (in  addition  to 
the  gasdynamic  Mach  number).  * 


l 


S!a 

»o“o 


**o  *ouo 


(1) 


For  an  interaction  region  of  length  L,  the  non- 
dimensional  numbers  are  defined  as 

L  L 

In/  l  dx  R  ■  /  r_  dx  <*> 

-u  «»  Jn  m 


When  R  >  >  1,  the  appropriate  measure  of  the 

intera&lon  is  the  parameter  S  defined  as 


S3  (Bo^o)/<Pui> 


where  the  spatial  average  (  )  is  over  the  elect¬ 
rically  conducting  portion  of  the  region  L.  For  a 
uniform  plasmold  of  length  a,  these  numbers  become 

I  •  '0B0‘/(>0U0-  Rm  =  VoV*  S  *  im¬ 
pulsed  magnetohydrodynamlc  flows  have  been 
examined  in  the  case  of  low  magnetic  ^Reynolds 
number  (R  <<  1)  and  weak  interaction  (I  <  I)1'*  . 
Because  of  the  low  interaction,  these  studies 
revealed  simple  current  flow  through  the  plasmold 
and  weak  magnetohydrodynamlc  deceleration. 

The  transverse  ionizing  shock  wave  which  forms 
the  front  of  the  plasmold  has  been  extensively 
studied  in  the  limit  of  infinitely  large  magnetic 
Reynolds  number1* ‘*T>.'  In  addition  to  the  exposition 
of  the  general  Ranklne-Hugonlot  conditions  for 
these  shocks*  it  has  also  been  demonstrated  tha 
such  shock  waves  can  be  reflected  as  well  as  trans¬ 
mitted  upon  encounter  with  an  externally  Imposed 
magnetic  field.  These  studies  also  showed  that  the 
electric  field  in  front  of  the  shock  must  be  self- 
consistently  determined  with  the  dynamical  state 
behind  the  shock  and  the  electrical  boundary 
conditions  imposed  upon  the  gas*. 

In  the  present  study,  we  examine  the 
magnetohydrodynamics  of  the  whole  plasmold  in  its 
encounter  with,  and  transit  through  an  externally 
imposed  magnetic  field.  We  show  that  under  con¬ 
ditions  of  strong  interaction,  hypervelocity 
plasmolds  can  possess  a  rich  variety  of  magneto- 
hydrodynamic  phenomena  including  magnetically 
reflected  shock  waves,  embedded  MHD  disconti¬ 
nuities,  and  significant  periods  of  transonic  flow 
within  the  plasmold.  In  particular,  we  reveal  the 
dynamics  of  reflected  and  transmitted  waves  through 
the  plasmold  in  both  the  low  and  high  magnetic 

Reynolds  number  regime.  We  reveal  the  behavior  of 
elcctrothermally  unstable  plasmolds.  We  show  that, 
in  general,  the  plasmold  is  not  delimited  by  the 
region  between  the  shock  front  and  contact  surface. 
Instead,  the  plasarald  develops  its  own  internal, 

evolving  structure  governed  by  the  mutual  inter¬ 

action  of  self-heating  and  Induced  self-fields. 

Shock-generated  hypervelocity  flows  of  this 
kind  are  subject  to  a  variety  of  nonideal 
phenomena.  These  Include  wall  interaction  effects 
(viscous  losses,  gas  leakage,  and  ohmic  voltage 

drops  in  boundary  layers),  thermal  radiation 
losses,  and  kinetle/lonlzatlon  relaxation  effects 
behind  shook  waves.  In  the  present  study  we  Ignore 
these  effects  and  examine  those  phenomena  which 
arise  specifically  from  the  magnetohydrodynamlc 


1 


interaction. 


of  the  plasaoid  is  than  given  by  Ohm's  law  and  the 
Maxwell  equation  in  the  MHD  approximation: 


In  Part  II  of  what  follows  we  present  a  mathe¬ 
matical  description  of  the  flow  of  quasi-one- 
diaensional  plasaoid a.  In  Part  III  we  examine  the 
dynamics  of  strong  interaction  plasaoids  with  an 
applied  aagnetic  field  but  at  low  aagnetic  Reynolds 
number.  In  Part  IV  we  siallarly  consider  strong 
interaction  plasatoids  but  at  large  magnetic  Rey¬ 
nolds  number.  In  Part  V  we  illustrate  the  behavior 
of  "transitional"  plasaoids.  These  are  flows  in 
which  the  plasaoid  enters  the  aagnetic  field  at 
relatively  low  values  of  interaction  paraaeter  and 
aagnetic  Reynolds  number.  As  a  result  of  self- 
Joule  heating,  however,  the  plasaoid  conductance  is 
elevated  as  it  progresses  through  the  field 
carrying  it  into  the  strong  interaction,  high 
aagnetic  Reynolds  number  regime. 

II.  Mathematical  Description 
A.  Conservation  Laws  for  Fluid  and  Electricity 


We  consider  a  quasl-one-diaenslonal  descrip¬ 
tion  of  the  gas  aoving  over  the  spatial  coordinate 
x  in  tiae  t.  Average  properties  (over  the  cross- 
section)  of  the  duct  describing  the  flow  are  its 

density  P ,  Internal  energy  C  ,  velocity  U,  and 

pressure  p.  The  magnetic  field  7  we  separate  into 
an  applied  field  B.  and  a  field  T.  Induced  by  the 
plaaaM  currents.  Trom  these  we  select  the  total 
mass,  momentum,  and  energy  densities 

£p,m  *  pO.  esptt  +  U2/*)  +  B^/2p0j 

as  the  state  W(x,t)  specifying  the  flow  at  any 

point  in  space  aind  time: 


W(x,t> 


The  fluid  conservation  laws  are 


<<*) 


J  =  -  UBQ)  (7) 


where  rj  a  is  the  magnetic  diffuslvity. 

External  Interaction  conditions  with  an 
extarnal  circuit  including  inductive  coupling  with 
the  applied  aagnetic  field  coil  are  required  to 
complete  the  description  of  actual  flow  situations. 
Rather  than  Include  such  circuit  detail  in  these 
illustrations  we  asstaie  that  the  external  circuit 
is  configured  so  that  an  electric  field  ?  * 
f(0,E,0)  is  maintained  within  the  interelectrode 
region  whose  magnitude  is  uniform  in  space  and 
given  by 


E  ■  «(UB> 


(9) 


where  a  is  a  "load"  paraaeter  (0  s«s').  For  a 
passive  external  circuit,  the  value  k  *  0 

corresponds  to  a  shorted  external  circuit;  for  k  i 
t  the  external  circuit  is  open  circuited. 

The  fluid  variables  P  ,p,T,  <r,U  are  nondlmen- 
sionallzed  by  the  values  Po*PO>To.  cq.Uq 
characterizing  the  interior  of  the  initial  plaawoid 
before  encounter  with  tlie_aagnetlc  field.  Mon 
dlaenslonal  space  and  tiae  x,'t  are  defined  in  teras 
of  x,  nondlaenslonallzed  by  the  plasaoid  length  a, 
and  t  by  a/U  q  .  The  nondiaensional  parameters 
governing  the  interaction  are  the  Mach  niaber  M, 
the  gas  heat  ratio  Y  ,  either  of  the  interaction 
parameters  I  or  S,  and  the  magnetic  Reynolds  nianber 
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aw  3F_ 

at  3x  +  — 


(5) 


In  the  above,  F  represents  the  convected  fluxes  of 
mass,  aoaentuaT  and  energy  while  E  contains  the 
Lorentz  force  and  power  associated  with  the  applied 
magnetic  field  and  the  Joule  dissipation.  Jhese 
are  expressed  In  tanas  of  current  density  J  and 
aagnetic  field  7Q: 


m 

0 

L  * 

m2/p  ♦  p  ♦  B^/24*0 

3* 

(Txb°)x 

m/p(e  ♦  p ♦  B2/2u0) 

J2/e-7M0xBo) 

In 

this  Illustration 

study 

we  assume  that 

kinetic  effects  are  confined  to  the  relaxation 
layer  at  the  shook,  end  further,  that  the 
relaxation  layer  is  thin  compared  to  the  overall 
thickness  of  the  plasaoid. 

For  the  sla^ple  geometry  (x,y,z)  of  Fig,  I,  the 
aagnetic  field  s  is  given  by  9lO,0,B),  the  electric 
field  by  J*  *  fio.E.O),  and  the  current  by  3“  « 
T(0,J.0).  The  description  for  the  near  fields  J.B^ 


B.  Boundary  and  Initial  Conditions 

In  the  limit  of  R  —  «o  ,  the  system  consisting 
of  Eqs.  (A)  and  (7) "is  fully  hyperbolic.  For 
finite- R  ,  the  system  Is  mixed  hyperbolic/parabolic 
with  embedded  regions  where  resistive  effects 
occur.  The  boundary  and  initial  conditions  which 
specify  the  interaction  problem  for  an  explosion 
generated  plasaoid  encounter  with  a  magnetic  field 
are  as  follows.  As  an  initial  condition  we  take  an 
Idealized  explosion  driven  flow  In  which  the 
plasaoid  of  given  breadth  a  occupies  the  hot  zone 
between  contact  surface  end  shock  front'  (Fig.  3). 
At  time  t  «  0  the  shock  front  Is  located  at  the 
edge  of  the  magnetic  field.  Over  the  tiae  scale 
for  the  dynaalcs  of  interest  the  shock  front  of  the 
plasaoid  will  run  continuously  into  the  quleseent 
driver  gas  while  the  backward  running  rarefaction 
continuously  runs  into  the  explosive  source.  Hence 
the  boundary  conditions  for  the  fluid  equations  are 
those  of  specified  explosion  and  quleseent  states 
at  the  boundaries  x  »  x  a  respectively. 

The  boundary  condition  for  the  induced  aagnetic 
field  B  froa  Eq.  (8)  Is  that  of  syaaetry  across 
the  ovefall  plasaoid  so  that  at  x  «  L  ,  x  a  -L 
which  lie  outside  the  region  of  any  current  flow  2 


BjC-Lj.t)  * 

The  applied  magnetic  field  B.  is  uniform  in  both 
space  and  time.  v 

C.  Numerical  Procedures 

The  solutions  to  the  initial-value  problems 
formulated  above  and  to  be  discussed  in  Parts  III, 
IV,  and  V  are  computationally  generated  with  second 
order  accurate  explicit  finite  difference  opera¬ 
tors.  The  hyperbolic  system  is  treated  with  the 
MacCormack  version  of  the  lax-Wendroff-Richtmyer 
operator10  .  For  the  space-time  grid  utilized, 
comparisons  were  made  with  the  analytically 
available  solutions  for  the  zero  and  infinite  R  , 
zero  interaction  limits.  At  the  extreme  pressuPe 
ratios  of  10  between  driver  and  driven  gas  for 
these  explosion  generated  plasmolds,  the  maximum 
variations  between  the  computationally  generated 
and  analytical  solutions  within  the  plasmold 
(expressed  as  a  fraction  of  the  analytical 
solution)  are  0.025  in  velocity,  0.04  in  pressure, 
and  0.08  in  temperature. 

Ill,  Interaction  at  Low  Magnetic  Reynolds  Humber 

We  now  proceed  to  the  first  of  several 
illustrations  of  the  foregoing  description.  We 
consider  first  the  strong  interaction  of  a  plasmold 
with  the  magnetic  field  but  at  low  magnetic 
Reynolds  number.  In  Fig.  2  the  kinematics  of  this 
situation  are  shown.  When  the  incident  plasmold 
encounters  the  magnetic  field,  the  leading  shock 
front  may  be  both  transmitted  and  reflected.  In 
the  case  of  reflected  fronts,  the  rear  (contact 
surface)  of  the  plasmold  subsequently  interacts 
with  the  reflected  shock  front.  This  colliding 
disturbance  then  radiates  a  fast  and  slow  reflected 
shock  (denoted  shock  II)  back  through  the  plasmold 
(which  consists  of  subsonic  flow  behind  the 
reflected  shock  front)  where  it  then  collides  with 
the  now  strongly  decelerated  shock  front. 

For  the  illustration  shown  here,  we  Select  a 
plasmold  with  Mach  number  M  x  1.64,  interaction 
parameters  I  x  20,  S  x  200,  and  Reynolds  number  R 
x  0.1,  Just  before  encountering  the  magnetic  field™ 
The  full  conditions  for  the  flow  are  given  in  Table 
I.  We  Impose  the  condition  that  the  electrical 
conductivity  is  spatially  uniform  within  the  high 
temperature  plasmold  (we  oonslder  electrical 
conductivity  functions  which  are  consistently 
coupled  to  the  gas  thermodynamic  state  in  Part  V). 
This  uniform  conductivity  distribution  is  achieved 
In  the  computations  with  the  model  conductivity 
function 
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reestablishes  high  velocity  flow  through  the 
plasmoid  and  reencounters  the  transmitted  shock 
front  which  has  been  decelerated.  During  this 
period  of  strong  wave  dynamics  the  current 
distribution  within  ‘the  plasmoid  is  strongly 
affected  (Fig.  5).  The  current  is  diminished  to 
very  small  values  during  the  period  of  plasmoid 
deceleration  behind  the  reflected  shock  I,  and  then 
returns  back  to  enhanced  levels  after  passage  of 
the  fast  reflected  shock  II.  The  low  magnetic 
Reynolds  number  of  the  plasmoid  allows  the  current 
to  diffuse  nearly  uniformly  throughout  its  breadth. 

IV.  Interaction  at  High  Magnetic  Reynolds  Humber 

We  next  consider  the  behavior  of  a  uniform  con¬ 
ductivity  plasmold  in  the  high  magnetic  Reynolds 
number  regime.  For  this  case  the  incident  plasmold 
has  an  interaction  parameter  I  x  20  as  in  the  pre¬ 
vious  illustration  but  a  magnetic  Reynolds  number 
R  x  5.  Correspondingly,  the  interaction  mxxber  S 
has  the  reduced  value  S  x  4.  The  full  conditions 
for  this  flow  are  given  in  Table  II.  The  encounter 
of  this  plasmoid  with  the  magnetic  field  is  similar 
to  that  of  Part  III.  The  features  peculiar  to  the 
higher  Reynolds  number  are  best  perceived  in  the 
current  distribution  of  Fig.  7  which  is  more 
nonuniform  compared  to  that  of  Fig.  8.  When  the 
current  levels  rise  behind  the  reflected  shock  II, 
they  do  so  by  directly  following  the  shock  until  it 
merges  with  the  transmitted  front.  The  current  max¬ 
imum  then  remains  at  the  shock  front  of  the  plas¬ 
moid.  As  a  result  of  decelerating  Lorentz  forces 
concentrated  inaedlately  behind  the  shock  front, 
the  shock  front  is  slowed  and  the  overall  breadth 
of  the  plasmoid  is  decreased  as  it  progresses 
through  the  magnetic  field.  This  is  in  contrast  to 
the  plasmold  dynamics  of  Parts  III  and  V. 

V.  Transitional  Plasmolds 

We  now  turn  to  consideration  of  plasmold  behavior 
with  a  coupled  electrical  conductivity  model.  In 
contrast  to  the  previous  illustrations  in  which  the 
conductivity  is  spatially  uniform  within  the  high 
temperature  plasmold  and  vanishes  outside,  we 
consider  a  conductivity  which  is  appropriately 
coupled  to  the  thermodynamic  state  of  the  gas.  As 
a  result,  local  regions  within  the  plasmold  can  be 
rendered  more  conductive  by  the  self-Joule  heating 
of  the  plasmold.  With  the  electrical  conductivity 
strongly  coupled  to  the  Joule  dissipation,  a 
plasmoid  in  the  low  I,  low  R^  range  can  evolve  into 
the  large  I,  large  R  range  as  it  progresses 
through  the  magnetic  find  and  experiences  further 
Joule  heating.  We  term  such  flows  "transitional* 
plasmolds. 

We  consider  a  conductivity  function  of  the  form 
r*‘  ■%n  +  <r«i  (10) 


which  effectively  switches  on  a  constant  conduc¬ 
tivity  within  the  plasmold  and  switches  the 
conductivity  off  outside  the  zone  in  which  the 
plasmoid  exists.  The  dynamics  of  the  interaetlon 
are  exhibited  in  Figs.  4,  5,  and  6.  When  the 
plasmoid  shock  front  encounters  the  magnetic  field, 
it  is  both  reflected  and  transmitted.  The 
reflected  shock  I  collides  with  the  eontaet  surface 
and  initiates  a  fast  reflected  shock  II  and  a  slow 
reflected  shock.  The  fast  reflected  shock  II 


In  the  above  ?  is  an  electrical  conductivity  of  a 
neutral  specieV  background  and  n  .  is  the  Coulomb 
conductivity.  We  use  as  a  summary  representation 
of  these  two  contributions  the  forms 


n 
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where  Z  ff  is  the  average  effective  ionic  charge 
and  r  .  T  ,  in  A,  n  >0  are  parameters  for  a  given 

8 as.  0  0 

The  conductivity  function  Eqs.  (10)-O2)  is 
dominated  by  the  partially  ionized  conductivity 
at  low  temperature  and  goes  over  to  the  CoulomB 

(fully-ionized)  conductivity  at  high  temperature. 
This  conductivity  function  has  the  property  Be/BTS 
0  over  the  entire  range  of  temperature.  Since 

there  are  no  thermal  energy  loss  mechanises  (which 
would  be  principally  radiative)  included  in  the 
model,  the  plasmold  is  unconditionally  electro- 
thermally  unstable11’11.  .  This  convective 
instability  is  simply  a  growth  of  temperature  non- 
uniformities  within  the  plasmold  due  to  intensify¬ 
ing  Joule  heating  resulting  from  growing  electrical 
conductivity. 

The  interaction  of  a  representative  transi¬ 
tional  plasmold  is  shown  in  Figs.  8-10.  This 

plasmold  has  interaction  parameters  1  *  1,  S  =  1 
and  magnetic  Reynolds  number  R  *  1  Just  before  it 
enters  the  magnetic  field.  The  complete  conditions 
for  this  flow  are  given  in  Table  III.  It  should  be 
noted  that  the  interaction  at  entry  into  the  mag¬ 
netic  field  is  considerably  smaller  than  the  inter¬ 
action  described  in  Parts  III  and  IV.  The  plasmold 
progresses  into  the  field  where  it  begins  to  self¬ 
heat  and  decelerate.  The  modest  interaction  at 
plasarald  entry  to  the  magnetic  field  does  not 
create  distinct  reflected  waves,  but  rather  a 
strong  and  continuous  deceleration.  Magnetic 
Reynolds  number  and  Interaction  parameter  grow 
slgnirioantly  has  the  plasma  is  heated.  At  time  ~ 

>  2.6,  the  magnetic  Reynolds  number  is  in  excess  of 
10  and  the  current  has  progressively  become  sheet¬ 
like.  It  should  be  noted  that  the  current  maxima 
no  longer  follow  lmsedlately  behind  the  transmitted 
shock  front.  Rather,  the  current  concentrates  in 
the  eleetrothermally  heated  zone  which  is  then 
convected  at  the  local  fluid  speed  rather  than 
radiated  at  the  shock  front  speed.  This 
decelerating  electrothermal  instability  is  then 
swept  up  by  collision  with  the  waves  initiated  by 
the  arrival  of  the  contact  surface  at  the  magnetic 
field  inlet.  A  feature  of  note  la  the  development 
of  reversed  current  flow  in  the  upstream  portion  of 
the  current  sheet  due  to  the  sharply  diminished 
magnetic  field  behind  the  current  sheet  at  large 

VI.  Summary  Remarks 

In  this  study  we  have  illustrated  significant 
purely  magnetogasdynamlc  phenomena  which  occur  when 
a  hypervelocity  pulse  of  plasma  ("plasmold”)  en¬ 
counters  an  applied  magnetic  field  under  strong 
interaction  conditions.  With  uniform  eleetrloal 
conductivity  within  the  plasmold  (and  vanishing 
electrical  conductivity  outside),  reflection  and 
transmission  of  the  plasmold  shook  front  are 
possible  coupled  with  strongly  nonuniform  current 
evolution  in  time.  Such  plasmolds  with  large 
magnetlo  Reynolds  numbers  have  current  distri¬ 
butions  (and  decelerating  Lorentz  forces)  concen¬ 
trated  lmsedlately  behind  the  shock  front.  As  a 
result  of  shook  front  deceleration,  these  plasmolds 
diminish  in  breadth  as  they  proceed  through  the 
magnetic  field.  kith  eleetrloal  conductivity 
within  the  plasmold  eoupled  to  its  thermodynamic 


state,  the  plasmoid  is  eleetrothermally  unstable 
and  creates  its  own,  evolving  region  of  enhanced 
electrical  conductance  which  carries  most  of  the 
current  and  is  convected  at  the  fluid  speed  within 
the  plasmoid.  The  shock  front  becomes  increasingly 
free  of  current  and  runs  progressively  farther 
ahead  of  the  unstable  current  structure  embedded  in 
the  plasmold  interior. 

Nonideal  phenomena  such  as  viscous  wall  layers, 
kinetie-relaxatlon  effects  behind  the  shock  front, 
and  thermal  radiation  losses  can  play  Important 
roles  in  these  strong  interaction  flows.  The  basic 
structure  of  the  magnetohydrodynamlc  interaction 
Itself,  however,  is  a  prerequisite  to  the 
description  and  understanding  of  these  additional 
modifying  effects. 
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Fig.  U.  Velocity  field  of  plasmoid  in  transit  through  an  applied  magnetic  field  under  strong  interaction 

(1  ■  30),  low  magnetic  Reynolds  number  (Ra  «  0.1)  conditions.  Shock  front  of  plasmoid  is  both  reflected 

from  magnetic  field  (a)  and  transmitted  into  magnetic  field  (b).  Reflected  wave  I  (a)  collides  with  con¬ 
tact  surface  (d)  and  initiates  a  fast  reflected  shock  II  (e)  and  slow  reflected  shock  (g).  Reflected  wave 

II  (e)  marges  with  transmitted  shock  front  (f).  Contact  surface  (h)  separates  conducting  from  nonconduct, 
lng  gas  behind.  Note  transitory  region  of  sharply  diminished  velocities  within  the  plasmoid  during  the 
period  of  reflection  and  rereflection  of  the  shock  front. 
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Mg.  5.  Current  density  distribution  for  conditions  of  Fig. 4.  Initial  peak  of  current  density  (a)  dim¬ 
inishes  sharply  as  plasmoid  velocities  diminish  behind  reflected  shock  front  1  (b).  After  passage  of 
reflected  shock  II  (e)  current  levels  rise  (d)  and  remain  uniform  through  plasmoid.  With  uniform  con¬ 
ductivity  model,  conductivity  is  uniform  between  shock  front  (e)  and  contact  surface  (f).  At  low  magnetic 
Reynolds  number  (RB  •  0.1),  current  is  diffused  nearly  uniformly  over  plaamoid. 
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Pig.  6.  Temperatu re  distribution  for  conditions  of  Fig.  U.  Othsrvlsa  uniform  temperature  of  plasnold  is 
affected  by  wave  dynamics,  and  most  significantly,  by  Joule  heating.  Weak  initial  current  levels  due  to 
reflected  shock  front  I  (a)  as  shown  in  Pig. U  lead  to  weak  elevations  in  temperature  due  to  Joule  heat¬ 
ing  (b).  Note  that  heated  region  within  plasnold  at  that  time  Is  convected  at  low  velocity.  Transmitted 
shock  front  is  sharply  decelerated.  Contact  surface  is  much  less  decelerated  during  this  period.  After 
passage  of  reflected  shock  II  (c) ,  Joule  heating  increases  significantly  thereby  elevating  plasmoid 
Interior  temperature  (b).  Following  passage  of  reflected  shock  II,  plasmoid  interior  is  convected  at 
velocities  ranging  between  that  of  front  and  that  of  contact  surface.  Trace  (e)  is  due  to  the  slow 
reflected  shock  wave  II. 


Pig.  7.  Current  distribution  for  a  strong  interaction  plasmoid  at  %  ■  5.  Dynamics  are  similar  to  those 
exhibited  in  Pigs.  U-6,  however  the  larger  Reynolds  number  permits  more  nonuniform  current  (and  induced 
magnetic  field)  distributions.  Current  pulse  initiated  by  encounter  of  plasmoid  shock  front  (a)  gives 
way  to  diminished  current  levels  la  low  velocity  region  behind  reflected  shook  front  I  (b).  Current 
levels  rise  and  shift  to  front  of  plasmoid  (e),(d)  after  passage  of  reflected  shook  front  II.  Note  that 
currant  concentrates  lamsdlatsly  behind  shook  front  In  contrast  to  that  of  Pig. 5  (whore  current  la  nearly 
uniform) and  that  of  Pig.  10 (where  currant  is  nearly  absent)  behind  shock  front. 
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Fig.  8.  Velocity  distribution  for  s  transitional  plasaoid.  Plasmid  enters  magnetic  field  with  1*1, 

Rn  ■  1>  3*1.  Modest  interaction  at  entry  to  aagnetic  field  does  not  generate  distinct  reflected  vaves ,  but 
rather  strong  and  continuous  deceleration.  Magnetic  Reynolds  number  and  interaction  parameters  grow 
during  the  course  of  transit.  Shock  front  proceeds  through  aagnetic  field  with  continuous  deceleration 
(a).  Once  the  Reynolds  number  reaches  values  Rg,  $  10,  the  current  has  become  sheet-like  with  a  deflag¬ 
ration-like,  intensifying  discontinuity  within  a  large  subsonic  tone  of  the  plasmid.  Contact  surface 
encounter  with  aagnetic  field  (b)  Initiates  downstream  running  waves  (c),  (d),  (e). 
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Pig.  9.  Current,  distribution  for  transitional  plasmid  of  Fig. 6.  Current  concentration  behind  shock  front 
at  entry  to  ascetic  field  grows  (a)  with  corresponding  heating  and  conductivity  enhancement  (reflected 
in  temperature  field  of  Fig. 10) .  Shock  front  propagates  ahead  (b),  but  becomes  increasingly  free  of 
current.  Electrotheraally  heated  tone  becomes  the  mat  highly  conducting  region  of  the  plasmid.  This 
heated  tone  is  convected  with  the  gas  and  strongly  decelerated.  Sharp  rise  in  temperature  and  current 
occurs  when  unstably  evolving  current  pulse  (which  is  convected)  collides  with  radiated  wave  (c).  Current 
pulse  is  then  swept  up  by  contact  surface  and  reaccelerated  (e) .  Rote  current  reversal  due  to  sharply 
reduced  aagnetic  field  in  upstream  portion  of  current  sheet  (f). 
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fig.  10.  Teapereture  field  corresponding  to  figs.  8,9  showing  developing  electrotherael  instability. 
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I.  INTRODUCTION 

In  this  paper  we  amine  sous  phenomena  of  interest 
in  megnetohydrodynaaic  (MHO)  chennel  flow  at  high 
■agnatic  Reynolds  number .  Such  flows  nay  be  created  by 
cheeical  explosions  and  the  ensuing  HNS  process  stay  be 
either  inherently  isisteady  1  or  quasi-steady.1 

In  the  unsteady  process  (Fig.  i>  the  flow  consists 
of  a  hot  plaaaui  ( “plaamold* )  formed  between  a  driven 
ionising  shock  wave  and  its  following  contact  surface. 
The  plaaaoid  Is  created  by  a  sudden  rales se  of  energy 
into  a  driver  section  which  is  in  contact  with  a  test 
gas  in  which  the  plaseold  propagates.  The  conducting 
plaseold  enters  a  region  in  which  an  externally  in posed 
■agnatic  field  Bl  and  electrodes  coupled  to  an 
asternal  circuit  exist  (Pig.  ').  The  plaana  conducts 
current  to  this  asternal  circuit  and  is  subject  to 
Lorents  forces  and  Joule  heating  as  it  propagates 
through  the  eegnetic  field.  If  the  explosion  drive  is 
a  chaeical  source,  such  a  pla^oid  will  be  of  the  order 
of  $-20  an  in  length  in  traversing  a  ■agnatic  field 
region  of  the  _<ydar  of  100  a  at  velocities  of  the 
order  of  10  as  The  plaaaoid  say  exist  at  pressures 
up  to  1  k  bar  and  energies  of  up  to  5  eV. 

In  a  steady  flow  device,  the  physical  schaaatie  is 
the  sane  as  in  Pig.  1  except  that  a  stagnation  chaaber 
exists  ahead  of  the  MB)  channel  in  tdiich  the  explosion 
driven  gases  are  processed.  This  configuration  yields 
a  hypersonic  channel  flow  which  is  qussi-steady  over 
the  flow  tine  scale  in  the  chennel. 

If  r  ,  P0,a0  are  the  characteristic  electrical 
conductivity,  ease  density  end  velocity,  the  flow  nay 
be  specified  by  an  interaction  nunber  i  and  aagnetle 
Reynolds  nunber  r  ( in  addition  to  the  gasdynanic  Mach 
nunber):  a 


In  Part  IV  we  consider  the  three-diaanaional  strong 
interaction  dynaaica  of  a  steady  hydronagnatic  channel 
flow  at  high  aaqnetic  Reynolds  n unbars .  The  quasi- 
one-diBensional ,  strong- interaction,  high  R  unsteady 
problea  la  treated  in  Rsf.  3.  * 

II.  Pnstaady  High  Magnetic  Reynolds  Nunber 


A.  Governing  Electrical 


Let  us  now  eonsider  the  fora  of  the  electrical 
equations.  The  Maxwell  equations  and  Ohn's  law  govern 
the  electric  field  I,  magnetic  field  •,  and  current 
density  7.  Let  t  he  the  axis  aligned  with  an  applied 
■agnatic  field  T°  and  let  x.y  bn  the  coordinates 
defining  the  plane  of  conduction.  The  current  vector 

then  becoaea  J  -  the  aagnetle  vector 

becoaes  *  -  »(0,0,I),  and  the  velocity  U  -  UIU^.U^.O) 

where  J  la  the  current  vector  noodlnons localised  on 
B°ypga.  »^0’  is  the  ■axlsna  value  of  the  applied 

aaqnetic  field.  U  is  the  velocity  nondlaanalonallsed 

on  the  duct  inlet  velocity  ug.  I  is  the  magnetic  field 
nondinepaionalised  on  the  max lmue  applied  aagnetle 
field  B  .  Ms  divide  the  aagnetle  UPld  into  the 
applied  "portion  7°  and  a  portion  7  1  which  is 
induced  by  the  currents  flowing  in  the  plaana.  The 
Maxwell  equations  and  Ohm's  law  then  aay  be  coabined  in 
the  induction  equation  fora  as 

ip  -  i  (U I  w".  •  Asv f;  &  ■  ££!)  • 

9t  dx  *  ay  Y  V9x*  9y2  ' 


M,o<rouo 


For  an  interaction  region  of  length  L,  the  nondiaen- 
» tonal  nuabers  are  defined  as 


is  /  i dx 
0 


/  r  dx 

-t  m 


When  s  >»  1,  the  appropriate  eeasure  of  the  Inter¬ 

action  %  the  parameter  f  defined  as 


.  ia  iSl!  *  ia  .  jL(5  s'01) -Au  af0>i  - 

3x  9x  3y  9y  J  9x  x  Y  Y 


In  Eq.  («),  n  ■  Vj/c  is  the  nondiaensional  eegnetic 
dlffusivlty,  Ra  is  the  aagnetle  Reynolds  nuabor  based 
upon  Ugb  ,  and  a,  and 


S  s  <b£/w0)  /  ,  pUZ) 


where  the  spatial  average  <  >  is  over  the  electrically 
oond acting  portion  of  the  region  L.  For  a  uni fore 

p  la  saw  id  of  loagth  a.  these  nuabers  bsceas 

I  *  V"«V  *■  •  V#V-  *  *  i,\- 

In  fart  IX  of  whet  follows  we  consider 
t  dimensional  electrical  conduction  <sidar  conditions 
of  high  R  for  the  unsteady  flow  of  a  plaaaoid  and  its 
encounter  *  with  the  applied  eegnetic  field  and  the 
electrode  tystan.  In  fart  III  we  consider  the 
two- dimensions!  steady,  hypersonic  flow  of  plasma  at 
high  aeqnotic  Reynolds  nuabor  through  a  constant  area 
reotaagelar  dust.  US  describe  the  currant  distribution 
is  the  plaMa  of  Farts  tx  and  III  under  the  conditions 
of  high  R  but  for  vanishing  inesractlon  I. 


The  current  density  components  are 

-!  dS11 

J  x  u-1  «■—  J  a  -u„ 

x  n>  y 

9y  1 


The  boundary  conditions  appropriate  to  Bq.  («>  are 
that  on  insulating  boundaries. 


while  on  conductors 


where  n 
aurfsca. 


is  the  nor 


coordinate  to  the  conductor 


T 


B.  The  Unsteady  flow  Ptawil 

In  thn  unsteady  flow  process,  a  plaaaoid  of  initial 
braadth  a  antara  the  eaqnetic  fiald  and  alaetroda 
ration  at  line  t  ■  0.  tha  front  of  tha  plaaaoid  moves 
at  tha  a hock- front  valocity  0  .  Tha  fluid  partlclaa 
bahind  tha  shock  front  and  beirend  tha  contact  aurfaca 
■ova  at  tha  plaaaoid  velocity  u  .  within  tha  plaaaoid 
(deaurkad  by  tha  shock-front  and  oontact  aurfaca),  tha 
electrical  conductivity  is  r  ,  while  outside  tha 
conductivity  has  tha  value  <r  «  cjj  . 


dovnetreaa  froa  the  duct  entry  froa  tha  driver,  it  tha 
station  lo  o  upatraaa  froa  tha  generator  inlet  ( x  • 
0.50  a)  tha  conditions  of  tha  flow  are 


Velocity  at  duct  centerline 
Teaperature  at  duct  centerline 
Mach  nuaber 

Boundary  layer  thickness 


9612  a/sec 
36,800  K 
3.  14 
5.5  aa 


at  tha  canter  of  the  generator  section  (x  •  0.60  a) 

those  values  are 


Tha  applied  ■agnatic  field  alonq  the  duet  axis  is 
assuaed  to  be  <ji van  by  the  following  aodal  profile 
function  s 

r 


velocity  at  duct  centerline 
Teaperature  at  duct  eentarline 
Mach  nua&er 

Bo undry  layer  thickness 


9520  a/sec 
36.930  K 
3.11 
6.5  aa 


ht  tha  exit  of  the  duct,  (x  ■  0.70  a)  these  values  are 


l+exp(-aL)-exp[-o(x+L/2)] 


(8) 


ll-exp(-«L/2)J2 

axp[a(x-L/2) ) 
(l-expt-oL/i)^ 


‘l<7 


where  B  is  the  aaxianaa  applied  field.  The  values 
selected1  are  (.  •  4h  and  a  •  2/h  where  h  is  tha  duet 
height. 


C.  Besaltlno  electrical  oistrlbutlona 

Given  an  external  load  condition,  and  a  piven 
plaaaoid  paadynaale  condition  we  seek  to  describe  the 
reaultinp  current  and  induced  aapnatic  field  distri¬ 
butions  as  the  plaaaoid  progresses  throuph  the  WO  duct 
and  panerator  section.  ue  shall  consider  conditions 
for  U_/U  -  4/3  ( oorrespondinp  to  stronp  shocks  in  noble 
gases!  and  initial  plaaaoid  breadths  ag/h  •  1.  The 
aaxisnai  nondiaenaiooal  current  flowinp  in  the  load 
circuit  is  ?  i  //(^  E?’">  where  I  is  tha  aaxianss  current 
par  isiit  length  flowing1  to  the  external  load. 

The  current  atreasUine  distributions  (and  corraapondinp 
induced  aapnatic  field  dlatribution I  for  a  aapnatic 
Reynolds  nuaber  ■  •  10  are  shown  in  rip.  2.  The  con- 

plate  dynaaics  of  interaction  between  the  panerator 
current  and  the  eddy-currents  Induced  by  the  aapnatic 
field  distribution  are  revealed  in  rip.  2. 


velocity  at  duct  centerline 
Teaperature  at  duet  centerline 
Mach  nuaber 

Boundary  layer  thickness 


9430  a/sec 
36,970  K 
3.07 
7.5  wm 


This  distribution  of  nonunifora  velocity  and 
teaperature  in  x  and  y  (with  corraapondinp  density  and 
conductivity  nonuni foralty)  reaultinp  froa  the  boundary 
layers  is  used  as  the  basis  for  the  electrical 
calculations  utilisinp  Bp.  (4).  The  nondiaensionalixed 
values  are  baaed  upon  thn  velocity  and  conductivity  at 
the  duet  inlet  (x  -  0).  the  Reynolds  nuaber  is  baaed 
upon  the  channel  height ,  h  and  has  the  value  R  •  7. 
The  velocity  and  teaperature  profiles  at  various 
stations  sloop  the  duct  are  shown  in  Tip.  4.  Since 
these  calculations  decouple  the  electricity  froa  the 
fluid  behavior  (weak  or  vanishing  interaction)  only  tha 
aapnatic  Reynolds  nuaber,  R  ,  is  a  relevant  electrical 
paraaetar.  " 


B.  Resulting  Blectrlcal  Distributions 

The  appropriate  electrical  equations  to  be  solved 
are  Eqa.  (4)  and  (5)  under  steady  stats  conditions. 
Boundary  conditions  ara  as  piven  in  Bqns.  (6)  and  (7). 
The  velocity  and  teaperature  profiles  (and  the 
corresponding  electrical  conductivity  profile)  will  be 
*  in  liar  to  those  shown  in  Pip.  3 1  however  we  shall 
reveal  the  affects  of  different  absolute  electrical 
conductlvllty  levels  (or  velocity  or  size  levels)  by 
varying  the  aapnatic  Reynolds  maber  for  the  various 
coaputed  results. 


III.  Steady  IWO  Channel  Plow  at  High  R 

h.  Channel.  Applied  Maenatlc  Plaid,  and  Plow 
Configuration 

we  now  consider  the  electrical  conduction  reaultinp 
froa  the  steady  flow  of  plasaa  at  high  aapnatic 
Reynolds  nuaber  in  a  constant  area  duct.  The  flowe 
exist  within  the  duct  gaenetry  and  applied  aapnatic 
field  distribution  shown  In  rip.  1  in  which  the 
panerator  electrode  length  is  equal  to  tha  duet  height. 
The  applied  aapnatic  field  is  piven  by  Bq..  (8). 

The  flow  field  has  been  calculated  as  that 
resulting  froa  a  reservoir  which  feeds  the  duet  with 
ergon  at  an  inlet  velocity  of  10  tae/eec.  an  internal 
energy  of  i?  HJ/kg,  a  pressure  of  10  k  bar,  and  a 
nealnal  electrical  conductivity  of  25,000  aho/e  at  the 
duct  inlet.  The  resulting  turbulent  boundary  layer 
interaction  is  considerable  as  tha  flow  proceeds  down 
the  duct.  These  flow  distributions  have  bean 
calculated  with  the  PTC  Research  Corporation  Q3D  faally 
of  oodesV 


In  Pip.  4  wa  exhibit  tha  open  circuit  and  loaded 
condition  for  a  aapnatic  Reynolds  nuaber  R  •  i.  it 
can  be  seen  that  there  ia  an  onset  of  ^ddy  cells 
induced  in  the  regions  of  ■agnatic  field  gradients  and 
a  convection  of  tha  current  pattern  downatreaa  of  the 
generator.  The  principal  generator  current,  however, 
is  still  confined  to  the  region  between  the  electrodes. 

In  Pig.  5  we  exhibit  the  sam  situation  but  for  a 
■agnatic  Reynolds  nuaber  R  •  7.  ns  note  that  consid¬ 
erable  current  flow  indi£ed  by  the  aiagnetlc  field 
gradients  exists  undar  open  circuit  conditions.  Mots 
how  the  upatraaa  eddy  cell  actually  couples  into  the 
generator  electrodes.  (Rider  load,  we  see  that  the 
generator  eurrent  exists  coooletelv  downstruao  froa  the 
electrodes. 

In  Pig.  6  the  conduction  field  for  a  aapnatic 
Reynolds  maber  Ra  *  (7.5  is  shown. 

again,  considerable  eddy  currant  activity  exists 
both  at  open  circuit  and  under  load.  The  generator 
eurrent  la  again  driven  out  of  tha  inters lectrode  pap 
and  exists  downstreaa  of  tha  electrodes. 


Tha  generator  electrodes  ara  located  at  0.6  a 


>>  i,  it  is  possible  to  show  thst  ths  transverse 
Lorentz  force  is  of  order  1/R  tines  the  axial  Lorentz 
force.  Hence,  cross  plane  aoaientum  is  weakly  generated 


xv.  Three-Dimensional  Strong  Interaction  I 
at  High  Magnetic  Reynolds  Number 

».  Flow  Structure  Prediction  in  Hydroeaqnetic  Channel 
flow 


In  short  burse  MHO  generator  channel  flow 
(tijaescale  100  u  s  -  1000  us)  utilizing  explosively 

generated  and  compressed  plasms  as  a  stagnation 
source,  the  flow  is  very  nearly  a  steady  channel  flow. 
At  the  plasma  conditions  of  electrical  conductivity, 
length  scale,  and  velocities  achievable  in  such  flows 
both  the  interaction  parameter  I  (ratio  of  Lorentz 
force  to  change  of  momentum)  and  magnetic  Reynolds 
number  R^  are  large.  Tor  ths  hydro* agnatic  flows 
typical  of  the  generators  characterized  by  the 
operating  conditions  of  Table  1  we  find  the  magnetic 
Reynolds  nirndter  R  to  be  about  20  and  the  interaction 
paraaMtar  X  to  be  about  IS /meter.  These  conditions 
indicate  that  the  deceleration  time  for  the  fluid  is  of 
the  order  of  the  diffusion  time  for  the  external 
magnetic  field  (consisting  of  applied  and  possible 
self -excited  combinations)  to  penetrate  through  the 
plaaaa.  One  may  therefore  expect  very  significant 
interaction  between  the  velocity  field  of  the  flow  and 
the  nonuniform  Lorents  forces  associated  with  the 
nomsiiform  magnetic  field  diffusion. 

In  Pig.  7  we  indicate  schematically  the  interaction 
of  the  magnetic  field  with  the  high  R  channel  flow. 
Awey  from  the  front  and  throughout  die  channel,  the 
lines  are  more  sharply  bent  through  the  magnetic 
boundary  layers  and  resemble  the  cross  plane  structure 
shown  in  Fig.  7.  At  high  magnetic  Reynolds  number  the 
magnetic  boundary  layers  are  thin  and  transverse 
gradients  in  y.z  are  much  stronger  than  axial  gradients 
in  x.  Correspondingly,  the  generator  current  is 
concentrated  in  the  magnetic  boundary  layers.  At  the 
channel  inlet  the  magnetic  field  is  most  strongly 
excluded  from  the  flovt  at  the  exit  the  field  has 
penetrated  most  deeply.  Ue  now  ask  how  a  hydraaagnetic 
flow  will  respond  in  its  velocity,  current,  Lorentz 
force,  and  power  generation  distributions  to  this 
fundamental  nonuniformity  situation  in  high  magnetic 
Reynolds  Humber,  strong  interaction  flow. 


TABLE  1 


Conditions  for  Illustrative  High  Magnetic 
Reynolds  Number  Channel  Plow 


Generator  Inlet  Conditions 
<U> 

T 

P 

a 

External  Magnetic  Field  ( »* ° 1 
Generator  Cross  Section  (h) 

R 

l" 


1.  The  Mathematical  Model 

The  fluid  mechanical  evolution  is  calculated  with 
the  STO  Research  Corporation  Q30  code  family 
described  in  *  and  *  .  In  the  version  of  the  Q30  family 
utilised  in  these  calculations  we  canter  attention  on 
the  axial  velocity  field  U  (x,y,s)  and  neglect  ths 
effects  of  the  cross  flow  ot  secondary  flow  velocity 
field.  Since  I  >>  t,  the  flow  is  dominated  by 
electrcaMgnatic  rather  than  fluid  mechanical  viscous 
effects,  although  turbulent  viscous  effects  are 
included  in  the  calculation.  Under  conditions  of  I  >> 


7  x  103  m/s 
2.0  x  104*K 
3  k.  bar 
S  x  10  mho/m 

)  10  T 
S  as 


by  Lorentz  forces  compared  to  axial  generation.  The 
primary  source  of  cross  plane  velocity  nonuniformities 
is  the  redistribution  of  axial  mass  flow  due  to  the 
nonuniform  deceleration. 

Let  us  now  consider  the  electrical  description. 
With  IT  -  (U  ,0,0),  the  three-dimensional  induction 
equation  can  be  expressed  as 

aff  _  3U  _  _ 

=  -B  x(B-7)Ux-7  xfiV  x  B)  (9) 


where  Xjis  the  unit  vector  in  the  x-direction  and  0  ■ 
(HqC)  is  the^  magnetic  diffusivity.  The  boundary 
conditions  for  7  are  as  follows.  As  can  be  seen  from 
Fig.  7  ,  specifications  of  the  boundary  data  for  a  on 
the  channel  perimeter  requires  solution  of  the  external 
(outside  the  channel)  as  well  as  internal  problem  for 
B.  As  an  alternative  to  solution  of  the  external 
problem,  we  approximate  the  boundary  conditions  for 
a  ,B  as 


w  / 


v. 


B  x 
z 


¥  +AB)exp(-Ax)  y  >  0 
±  (Bg +AB)exp(- Ax)  y<0 

[  1  *  exp(-Ax)  ] 


CO) 


on  y  -  *7 


B  a  0 

y 

B^  =  B^j  +  AB  exp(-Vx) 


(i  n  2 

where  A  -  /  [j \  ^2  is  the  decay  term  for  the  fundamental 
diffusion  mode  in  a  square  duct,  8  is  the  far-field 
external  value  of  the  magnetic  field  (totally  in  the  z 
direction),  and  A  a  is  the  maximum  compression  of  the 
magnetic  field  at  the  generator  inlet.  It  can  be  seen 
that  the  approximate  model  (Eq.  (10)  possesses  the 
correct  limiting  values  at  x  —  oo .  For  R  >>  1,  the 
values  at  x  —  0  are  also  close  to  the"  full-field 
exclusion  solutions.  The  modal  is  least  accurate  for 

conditions  of  R  «  1. 

m 


For  a  the  boundary  conditions  are  B  •  z  p.  X/2  on 
z  •  *  w/J  where  I  is  the  total  current *per  unit  depth 
flowing  in  the  generator  in  the  transverse  direction. 


C.  Magnetic  Diffusion  and  Flow  Interaction 

ue  now  consider  the  predictions  for  flow  in  the 
generator  channel  of  Fig.  7  under  the  conditions  of 
Table  1.  For  these  conditions  we  set  the  generator 
current  at  the  nominal  value 


Z  -  (1-K)  w(cu  B^) 

with  K  •  1/2  and  v  •  the  channel  width.  Ue  assume  that 
no  field  has  penetrated  the  plasma  at  the  electrode 
edge  x  •  0.  (In  general,  one  will  expect  a  starting 
magnetic  boundary  layer  thickness  at  the  channel 
entrance  dependent  upon  the  plasma  nozzle  geometry  and 
the  fringe  field  location.) 

The  distribution  of  the  transverse  current  2  and 
the  magnetic  field  component  a  for  various  ¥ixlel 
stations  in  the  channel  are  shown*ln  Fig.  8  .  (There 
is  an  axially  directed  eddy  current  J  associated  with 
the  fields  8^  ar.i  Sy  which  is  not  shoA>  here.'  it  x  • 


1.4  cm  the  magnetic  boundary  layers  are  exceedingly 
thin  Is  3  on).  They  have  grown  to  1.4  ca  it  x  »  14  ct> 
and  2.9  cm  at  x  *  28  cm.  The  Inward  diffusing  magnetic 
field  component  B  i*  initially  concentrated  along  the 
walla  parallel  to  the  external  magnetic  field  (elec¬ 
trode  walls;  see  Pig.  1),  but  the  z  component  on  the 
side  walls  builds  up  as  the  flow  proceeds  down  the 
channel.  The  transverse  current  J  it  initially 
concentrated  along  the  sidewalls  and  diffuses  Inward 
from  these  walls. 

In  Fig.  9  we  observe  the  impact  of  these 
distributions  of  current  and  magnetic  field  on  Lorentx 
force  and  axial  velocity.  In  the  early  stages  of  the 
flow  evolution  the  maximtas  axial  Lorentz  force  is  an 
the  corners.  This  is  because  the  maximtas  field  B 
exists  on  the  electrode  walls  parallel  to  the  z  axis 
and  maximtas  current  J  exists  on  the  sidewalls  parallel 
to  the  y  axis.  The  Maximtas  product  of  the  two  there¬ 
fore  exists  in  the  corners.  Throughout  the  entire 
evolution,  the  maximtas  lorentz  force  is  concentrated  in 
the  magnetic  boundary  layers  along  the  sidewalls. 

The  Impact  on  the  velocity  distribution  is  strong 
and  significant.  The  initially  flat  velocity  profile 
is  rapidly  decelerated  in  the  sidewall  layers.  By  12 
cm  the  flow  is  actually  reversed  in  the  corners  and  the 
regions  of  reverse  flow  spread.  Because  of  the 
deceleration  of  the  wall  regions  by  the  inward 
diffusing  Lorentx  force  the  core  flow  is  accelerated. 
Ftom  an  initial  value  of  7  km/s  the  core  flow  is 
accelerated  to  nearly  10  km/s  by  x  •  28  as.  Since 
there  is  little  Lorentx  force  in  the  core  region,  this 
high  velocity  region  is  not  effectively  coupled  by  the 
generator . 

It  should  be  noted  that  once  reverse  flow  occurs 
(stalled  HKD  boundary  layers)  the  present  prediction  of 
the  actual  velocity  field  in  the  wall  layer  region 
cannot  be  considered  correct  since  upstream  feedback  is 
prohibited  by  the  method  of  calculation.  The  present 
calculation  is  only  definitive  in  its  prediction  that 
stall  will  indeed  occur.  In  actuality  we  conjecture 
that  the  flow  under  these  conditions  will  appear  as 
shown  in  Fig.  10  in  which  recirculation  eddies  are 
trapped  within  the  boundary  layer  region  of  intense 
Lorentz  force  while  the  core  flow  jets  through  outside 
the  region  of  the  eddies.  In  such  a  stalled  flow, 
there  is  little  Lorentz  force  coupling  to  the  high 
maaentvn  jet  in  the  central  region  of  the  channel. 

As  a  demonstration  that  stall  is  not  a  peculiarity 
of  the  corner  region,  we  show  the  Lorentz  force  and 
velocity  profiles  at  t*  21  m  for  a  higher  interaction 
aituation  (  <t  -  2  x  10*  mho/m)  but  otherwise  conditions 
identical  to  Table  1  (Fig.  11).  In  Fig.  12  we  see  that 
the  entire  sidewall  has  reversed  flow  under  these 
conditions. 

In  Fig.  13^  we_  show  the  distribution  of  power 
production  •(  J  *  E)  over  the  cross  section,  power  is 
only  produced  in  the  mid-wall  layer  regions  and  is 
actually  consumed  in  the  corner  and  sidewall  layers 
where  reversed  flow  takes  place.  The  overall  power 
extraction  per  unit  length,  dP/dx,  is  defined  as  the 
integral  of  the  power  density  over  the  channel  cross 
section » 

»  -  /'  ft 7  •  E)»i2x 

dx  A 

The  ideal  power  attraction  for  uniform  flow  is 


(IE)  =  K(l-K)<ru2B2A 

'  dx 'ideal  0  5 


channel  IS  at  long,  only  10%  of  the  ideal  power  is 
available  and  for  a  28  cm  channel  only  27%  is 
available.  These  results  must  be  tempered  by  the 
observation  that  beyond  12  ®t  the  boundary  layers  have 
begun  to  stall;  it  is  not  likely  that  power  extraction 
will  be  enhanced  by  realistic  inclusion  of  the  recir¬ 
culation  zones  under  stalled  conditions.  Hence,  the 
power  extraction  levels  beyond  x  •  12  cm  will  likely 
not  be  reached. 
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In  Fig.  14  we  exhibit  the  ratio  of  tr.ase  quantities  for 
the  conditions  of  Table  i.  It  can  be  seen  that  for  * 
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Fig.  1.  Schema? i  :  of  explosion  driven  plesmoid  flow. 
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Flq.  2. 
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Indue*!  aagnetic  field  iaolavelt  (currant  danalty  streamlines)  at  R  -  10  in  unsteady  passage  of 
explosion-driven  plasaoid.  At  tins  t  «  0,  tha  shock  front  Is  at  the  "edge"  of  tha  applied  aagnetic  field. 
Upper  figure  in  each  tine  sequence  indicates  plasaoid  position.  Lower  figure  indicates  aagnetic  field 
lsolevels.  Note  intense  eddy-current  calls  in  (a)  and  (a)  when  plasaoid  is  in  region  of  applied  aagnetic 
field  gradient  and  out  of  range  of  electrodaa.  tn  (bl,  tha  generator  current  is  concentrated  in  a  thin 
tone  behind  the  shock  front.  In  'cl.  the  plasaoid  has  progressed  downstream  to  that  the  generator  current 
is  strongly  converted  downstream.  Ir.  'd).  the  generator  current  and  iovnstream  eddy  cell  interact 
together . 
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Pig.  3.  Velocity  and  teaperature  profile*  between 

electrode*  in  turbulent  «u parsonic  plaana  flow 
in  an  MHO  generator  duct. 
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Pig.  S.  Induced  aagnetic  field  isolavels  I*11  for 
steady  nagnetohydrodynaaic  flow  with  the 
profile*  of  Pig.  3  and  •  J.,  In  (a),  the 
generator  is  at  open  circuit  (_/•  01  end  in  (b) 
the  generator  la  under  load  ( I  •  1 ) . 
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Pig.  4.  Induced  aangetic  field  iaolevel*  B* 1  for 

•teady  nagnetohydrodynaaic  flow  with  profile* 
•taller  to  those  of  rig.  3  and  R^  •„  1.  In  (a), 
the  generator  1*  at  open  circuit  «  0)  and  in 
<bl  the  generator  i*  under  load  (/•  1). 
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Pig.  6.  Induced  aagnetic  field  iaolevel*  B*1*  for 

steady  nagnetohydrodynaaic  flow  with  profiles 
siailar  to  those  of  Pig.  3  and  R  *  17. 5*  In 
(a),  the  generator  is  at  open  circuit  (/^  0) 
and  in  (b)  the  generator  is  under  load  (  I  -  i) 
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Pig.  7.  Magnetic  field  nonuntfomitles  in  nigh  Magnetic 
Reynolds  y  vaster  channel  flow. 
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rtf.  12.  Velocity  distribution  at  x  *  21  cm  correspond' 
inp  to  s  niftier  Interaction  ease  ( <r  •  2  s  10* 
■no/a I  but  otherwise  conditions  as  in  Table  i. 
Note  that  entire  sidewall  flow  is  reversed  in 
this  case. 
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Ftp.  14.  Ratio  of  power  extraction  per  unit  lenptb  to 
ideal  power  extraction,  wall  layers  stalled 
over  nost  of  channel  lenpth. 


